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ABSTRACT 

Previous observations of the Ap star HD 32633 indicated that its magnetic field was unusually 
complex in nature and could not be characterised by a simple dipolar structure. Here we 
derive magnetic field maps and chemical abundance distributions for this star using full Stokes 
vector (Stokes IQUV) high-resolution observations obtained with the ESPaDOnS and Narval 
spectropolarimeters. Our maps, produced using the InversIO magnetic Doppler imaging 
(MDI) code, show that HD 32633 has a strong magnetic field which features two large regions 
of opposite polarity but deviates significantly from a pure dipole field. We use a spherical 
harmonic expansion to characterise the magnetic field and find that the harmonic energy is 
predominately in the £ = 1 and i = 2 poloidal modes with a small toroidal component. At the 
same time, we demonstrate that the observed Stokes parameter profiles of HD 32633 cannot 
be fully described by either a dipolar or dipolar plus quadrupolar field geometry. We compare 
the magnetic field topology of HD 32633 with other early-type stars for which MDI analyses 
have been performed, supporting a trend of increasing field complexity with stellar mass. 
We then compare the magnetic field topology of HD 32633 with derived chemical abundance 
maps for the elements Mg, Si, Ti, Cr, Fe, Ni and Nd. We find that the iron-peak elements show 
similar distributions, but we are unable to find a clear correlation between the location of local 
chemical enhancements or depletions and the magnetic field structure. 
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1 INTRODUCTION 

The magnetic chemically peculiar (Ap/Bp) intermediate-mass stars 
are main sequence objects of spectral types A and B exhibiting a 
number of unusual properties when compared to normal stars of 
similar spectral types. Their defining characteristic is the presence 
of strong, globally organised surface magnetic fields with a strength 
ranging from a few hundred G to a few tens of kG. In the first ap¬ 
proximation these fields have a geometry of a dipole inclined with 
respect to the stellar rotational axis (Donati & Landstreet 2009). 
The presence of such strong magnetic fields leads to a grossly non¬ 
solar atmospheric chemistry and to the formation of high-contrast 
horizontal (spots) and vertical (stratification) chemical abundance 
inhomogeneities. As a consequence of this non-uniform surface 

* Based on observations obtained at the Canada-France-Hawaii Telescope 
(CFHT) which is operated by the National Research Council of Canada, 
the Institut National des Sciences de l’Univers of the Centre National de la 
Recherche Scientifique of France, and the University of Hawaii. Also based 
on observations obtained at the Bernard Fyot Telescope (TBF, Pic du Midi, 
France) of the Midi-Pyrenees Observatory, which is operated by the Institut 
National des Sciences de l’Univers of the Centre National de la Recherche 
Scientifique of France. 


structure, Ap/Bp stars exhibit a profound spectroscopic and photo¬ 
metric rotational variability, even though the intrinsic configuration 
of spots and magnetic field remains constant on the time scale of 
decades. 

Detailed empirical information on the geometry of magnetic 
and chemical surface structures in Ap/Bp stars is essential for un¬ 
derstanding the physical mechanisms governing evolution of the 
global magnetic fields in stellar interiors (Braithwaite & Nordlund 
2006; Duez & Mathis 2010) and for testing theories of chemical 
segregation by radiative diffusion (LeBlanc et al. 2009; Alecian 
& Stift 2010). Aiming to collect this information for a meaning¬ 
ful sample of Ap/Bp stars, we are reconstructing magnetic and 
chemical maps of their surfaces from high-resolution spectropo- 
larimetric observations (Kochukhov et al. 2004, 2011, 2014, 2015; 
Kochukhov & Wade 2010; Silvester et al. 2014a; Rusomarov et al. 
2015), with a special emphasis on the interpretation of line pro¬ 
file variability in all four Stokes parameters. In the this paper we 
present the first such analysis for the Ap star HD 32633. 

The star HD 32633 (HZ Aur, HIP 23733) has an effective tem¬ 
perature of 12,000-13,000 K and is of the peculiar subclass SiCr 
(Renson & Manfroid 2009). Early observations by Babcock (1958) 
measured the magnetic field as a function of phase; this author 
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noted that the magnetic field was strong, reversing in sign approx¬ 
imately once a week and that the field behaved erratically at the 
negative polarity. Later observations by Borra & Landstreet (1980) 
and Leone et al. (2000) found that the longitudinal field variation 
as a function of rotational phase could not be fit with a simple si¬ 
nusoidal curve. Instead, the curve appeared non-harmonic, with ex¬ 
trema separated by 0.6 of the rotational period. These early obser¬ 
vations indicated that the magnetic field of HD 32633 departs from 
a simple dipolar structure. It was later suggested by Glagolevskij & 
Gerth (2008) that the magnetic field of HD 32633 could be modeled 
assuming two dipoles at opposite regions of the star, located near 
the rotational equator. HD 32633 does not show signs of pulsational 
variability in either the roAp (Joshi et al. 2006) or SPB (Aerts et al. 
2014) frequency ranges. 

These observations establish HD 32633 as an interesting target 
for mapping using magnetic Doppler imaging (MDI), from which 
we can obtain detailed information on the magnetic field geometry 
and thus better understand the unusual non-sinusoidal longitudinal 
field phase curve behaviour of this star. Furthermore, the pioneering 
MuSiCoS four Stokes parameter observations of Ap stars by Wade 
et al. (2000) demonstrated that HD 32633 shows the strongest lin¬ 
ear polarisation signatures among the Ap/Bp stars included in their 
survey. More recently, Silvester et al. (2012) collected a superb new 
data set of full Stokes vector observations of several Ap/Bp stars, 
including HD 32633, using the ESPaDOnS and Narval instruments. 

Spectropolarimetric data in four Stokes parameters can be suc¬ 
cessfully interpreted with the magnetic Doppler imaging method, 
in particular with the INVERS10 code (see Piskunov & Kochukhov 
2002; Kochukhov & Piskunov 2002) optimised for simultaneous 
reconstruction of the vector magnetic maps and chemical abun¬ 
dance distributions. A remarkable result of the application of MDI 
to four Stokes parameter spectra of Ap stars was the discovery that 
the magnetic topologies of at least some Ap stars depart signifi¬ 
cantly from low-order multipolar configurations and contain impor¬ 
tant small-scale components (Kochukhov et al. 2004; Kochukhov 
& Wade 2010; Silvester et al. 2014a). However, after a 2 CVn 
(HD 112413) studied by Silvester et al. (2014a), HD 32633 is only 
the second Ap star to be mapped using the modern high-resolution 
Stokes IQUV data of Silvester et al. (2012) and as such provides 
an exciting further insight into how the magnetic field topologies 
of Ap stars might deviate from the canonical dipolar geometry. In 
addition to investigating the magnetic field topology we also take 
advantage of the time-resolved spectropolarimetric data to map the 
chemical abundance structures on the surface of HD 32633 in order 
to examine correlations between the magnetic field topology and 
the locations of any abundance enhancements or depletions. 

The paper is organised as follows. Section 2 briefly describes 
the observational data. Section 3 discusses derivation of stellar pa¬ 
rameters. Section 4 summarises the procedure of magnetic Doppler 
imaging. The resulting magnetic and chemical abundance maps are 
presented in Sections 5 and 6 , respectively, and are discussed in 
the context of other MDI studies and theoretical investigations in 
Section 7. 


2 SPECTROPOLARIMETRIC OBSERVATIONS 

Observations of HD 32633 were obtained between 2006 and 2010, 
with both ESPaDOnS and Narval spectropolarimeters in all four 
Stokes parameters (Stokes IQUV). The full details of the observa¬ 
tions, including the observing log, are reported by Silvester et al. 
(2012). Briefly, the spectra have a resolving power of 65,000, cover 


the wavelength range 369-1048 nm, and have a typical signal-to- 
noise ratio of 600:1. A total of 20 full Stokes IQUV observations 
are available for HD 32633, with an additional Stokes IV obser¬ 
vation. Observations sample on average every 0.1 of a phase, with 
closer sampling at certain phases. As described by Silvester et al. 
(2012), the reduction of the Stokes IQUV spectra was performed 
on-site at the respective observatories using the Lib re-Esprit re¬ 
duction package (Donati et al. 1997). Normalisation of the spectra 
was performed order by order using an IDL code optimised for 
Ap stars. The rotational phases were computed according to the 
ephemeris of Adelman (1997). 


3 FUNDAMENTAL STELLAR PARAMETERS 

We obtained an initial estimate of the atmospheric parameters of 
HD 32633 by applying the calibrations by Netopil et al. (2008) to 
the uvbyf3 (Hauck & Mermilliod 1998) and Geneva (Rufener 1976) 
photometric observations of this star. Both sets of photometric ob¬ 
servations indicate T e ff = 12500-13000 K and log# = 4.0-4.3. 
Guided by these estimates, we computed a model atmosphere struc¬ 
ture for T e ff = 13000 K, log g — 4.0 and [M/H] = +1.0 dex us¬ 
ing the LLmodels code (Shulyak et al. 2004). Using this model at¬ 
mosphere we then determined abundances for about a dozen chem¬ 
ical species, including He, Mg, Si, Fe-peak and several rare-earth 
elements, by fitting the Synmast (Kochukhov et al. 2010) theoret¬ 
ical spectra to the phase-averaged Stokes I spectrum of HD 32633. 
The resulting individual chemical abundances were adopted for cal¬ 
culation of a new LLmodels atmospheric grid covering the rele¬ 
vant T e ff and log g ranges. In these calculations we took into ac¬ 
count modification of the line opacity due to Zeeman splitting and 
polarised radiative transfer (Khan & Shulyak 2006) assuming the 
field modulus of 10 kG (Glagolevskij & Gerth 2008). 

The final effective temperature was determined by fitting the 
spectral energy distribution (SED) in the optical (Geneva photome¬ 
try converted to absolute fluxes and spectrophotometry by Adelman 
et al. (1989)) and near-IR (2MASS fluxes) wavelength regions, us¬ 
ing the LLmodels code to model the SED. HD 32633 is known 
to be somewhat reddened, with E(B — V) ranging from 0.05 to 
0.12 according to different estimates the literature (e.g Adelman & 
White 1980). Here we adopted E(B — V) — 0.10 =L 0.02 as a 
compromise value. Figure 1 illustrates comparison of the observed 
SED with the model fluxes corresponding to T e ff = 12800 K and 
log# as 4.15. The latter value of the surface gravity was estab¬ 
lished by comparing Synmast calculations to the observed pro¬ 
files of hydrogen Baimer lines. As evident from Fig. 2, our choice 
of atmospheric parameters yields an excellent fit to the Ha, H /3, 
and H 7 lines in the spectrum of HD 32633. 

As a by-product of the SED fitting we determined an angular 
diameter and converted it to the stellar radius, R = 1.99 ±0.2572©, 
using the parallax tv = 5.57 d= 0.67 mas (van Leeuwen 2007). Our 
radius value agrees well with previous estimates by Kochukhov & 
Bagnulo (2006, 1.92 d= 0.3272©) and Leone et al. (2000, 2.472©). 

A summary of fundamental stellar parameters determined or 
adopted in our study is given in Table 1 . 


4 MAGNETIC AND CHEMICAL INVERSIONS 

An application of four Stokes parameter MDI requires knowledge 
of several other quantities in addition to the fundamental stellar pa¬ 
rameters determined above. In particular, we need to known the 
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UV and optical 
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Figure 1. Comparison between theoretical spectral energy distributions (lines) and observations (symbols) in the UV, optical and near-IR spectral regions. The 
two theoretical SEDs correspond to the calculations for T e ff = 12, 800, log# = 4.15 including effects of individual non-solar abundances of HD 32633 and 
magnetic held (solid line) and standard model (solar abundances, no magnetic held) for the same T e ^ and log g (dashed line). 
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Figure 2. Comparison between the observed (dark curve) and computed 
(bright curve) prohles of the Ha, H (3, and H7 lines in the spectrum of 
HD 32633. 


Table 1. Fundamental parameters of HD 32633. 


Parameter 

Value 

Reference 

7"eff 

12800 ± 500 K 

this study 

log g 

4.15 ±0.1 

this study 

R 

1.99 ±0.25 Rq 

this study 

v e sin i 

18.0 ± 1.0 kms- 1 

this study 

i 

80° ± 5° 

this study 

© 

90° ± 5° 

this study 

Prot 

6.43000 d 

Adelman (1997) 


projected rotational velocity v e sin z and to specify the inclination 
(tilt with respect to the line of sight) and azimuth (orientation in 
the plane of the sky) angles of the stellar rotational axis. A review 
of the literature did not reveal a definitive value for the inclination 
angle of HD 32633. Therefore to constrain the orientation of the 
rotational axis we ran a grid of inversions by varying both the incli¬ 
nation and azimuth angles. We determined a suitable value for both 
angles by Ending the pair of angles which resulted in the small¬ 
est deviation between the model spectra and the observational data. 
Similarly, to determine an optimal v e sin z, we started with the pre¬ 
viously determined value of 19.0 =b 2.0 kms -1 (Silvester et al. 
2012 ), performed inversions for a range of v e sinz values and se¬ 
lected the value which provided the smallest deviation between the 
observations and model profiles. 

Through this process we found that the following parameters 
resulted in the best fit between the model and observations: inclina¬ 
tion i m 80° =b 5°, azimuth angle of 90° zb 5° and v e sin i = 18.0 zb 
1.0 km s- 1 . We note that our inclination angle agrees with the value 
of this parameter estimated by Borra & Landstreet (1980) (77°) and 
by Glagolevskij & Gerth (2008) (84° zb 5°). The rotational period 
from Adelman (1997) together with the stellar radius and v e sinz 
determined in our paper imply an inclination angle close to 90°. 
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Figure 3. Surface magnetic field distribution of HD 32633 derived using Stokes IQUV computed with InversIO using a spherical harmonic description of 
the field topology. The spherical plots (at inclination angle of i = 80°), show distributions of: a) the field modulus, b) the horizontal field, c) the radial field 
and d) the field orientation. Each column corresponds to a different phase of rotation (0.0, 0.2, 0.4, 0.6 and 0.8) and the colour bars on the right side indicate 
the magnetic field strength in kG. 


The methodology of magnetic Doppler imaging employed in 
this study is similar to that used in the mapping of the Ap a 2 CVn 
(Silvester et al. 2014a, b). We rely on the Invers 10 code (Piskunov 
& Kochukhov 2002) to obtain best-fitting magnetic and chemical 
distributions. A regularisation function is used to facilitate conver¬ 
gence to a global chi-square minimum and guarantee that we find 
the simplest surface distribution that fits the data. A new feature of 
the magnetic field reconstruction is that we parameterise the mag¬ 
netic field topology in terms of a general spherical harmonic expan¬ 
sion (up to the angular degree i= 10) as described by Kochukhov 
et al. (2014). However, for comparison with our previous studies, 
we also performed inversions with discrete pixel maps for each of 
the three magnetic field components. In both cases the choice of 
regularisation parameter is essential. For the maps presented in this 
paper, the value of regularisation was chosen so that it gave the 
lowest total discrepancy between observations and the model and 
yet still reproduced the Stokes profiles without fitting to a signif¬ 
icant amount of noise. An additional constraint was chosen, such 
that the total regularisation must be at the minimum only a factor 
of 10 times smaller than the total discrepancy. After this limit the 
improvement to the discrepancy becomes increasingly smaller, but 
the map will continue to become increasingly patchy. For most of 


the inversions, a stepwise approach was taken with the code us¬ 
ing five steps of regularisation (decreasing in value by a factor of 
three) during the course of one inversion run. The result is five sets 
of fits between observations and the model and five resulting maps 
(each corresponding to a different regularisation parameter). This 
allowed us to determine the most suitable value of regularisation to 
use for a given inversion run. 

Many metal lines in the ESPaDOnS and Narval spectra of 
HD 32633 exhibit polarisation signatures suitable for modelling in 
four Stokes parameters. Previously we found that simultaneous in¬ 
terpretation of line profile variability of several chemical elements 
with diverse surface distributions yields the most robust magnetic 
maps (Silvester et al. 2014b). Taking these results into account, we 
selected 10 mostly unblended lines of Sill, Crll, Fell and Ndm 
for the simultaneous magnetic and abundance mapping in Stokes 
IQUV parameters. In addition, 6 lines of Mgll, Sill, Till and 
Ni II with comparatively weak Stokes QU signatures were used for 
reconstruction of the surface abundance distribution of these ele¬ 
ments from the Stokes IV spectra adopting a fixed magnetic field 
geometry. The atomic parameters of 16 spectral features employed 
in the MDI of HD 32633 are given in Table 2. For each spectral 
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Table 2. Atomic lines used for the mapping of HD 32633. The log gf val¬ 
ues are those provided by the Vienna Atomic Line Database (VALD Kupka 
et al. 1999). 


Ion 

Wavelength 

(A) 

log af 

Magnetic and abundance mapping 

Si 11 

6347.109 

0.170 

Cr 11 

4558.650 

-0.449 


4558.650 

-4.950 


4588.199 

-0.627 


4824.127 

-0.970 

Fe 11 

4583.829 

-1.860 


4923.921 

-1.320 


5018.436 

-1.220 


5169.028 

-1.250 

Nd in 

5050.695 

-1.060 


5102.428 

-0.620 

Abundance mapping 

Mg 11 

4481.126 

0.740 


4481.150 

-0.560 


4481.325 

0.590 

Si 11 

4130.872 

-0.824 


4130.894 

0.563 


5055.984 

0.530 

Tin 

4290.215 

-0.870 


4563.757 

-0.690 

Ni 11 

4067.031 

-1.834 


interval we also included relevant minor blends, which are omitted 
in this table. 


5 MAGNETIC FIELD STRUCTURE 

To characterise the magnetic field topology of HD 32633, maps 
were derived using the spherical harmonic expansion version of 
InversIO (Kochukhov et al. 2014) and, for an addition test, us¬ 
ing the discrete surface element version of the code (Piskunov & 
Kochukhov 2002). The final map discussed below was obtained 
using the spherical harmonic parameterisation. In this inversion we 
obtained the best fit between the observed Stokes IQUV parame¬ 
ters and model profiles by allowing the code to fit harmonics up to 
£ = 10 , including a toroidal component and allowing independent 
radial and horizontal poloidal field components. In other words, the 
three sets of spherical harmonic coefficients a^, m , and 7 
specifying the surface field geometry (see Kochukhov et al. 2014) 
were allowed to vary independently. The resulting magnetic field 
structure is presented in Fig. 3 and the corresponding fits to obser¬ 
vations are illustrated in Figs. 6 and 7. 

Our MDI modelling indicates that the field structure of 
HD 32633 is dominated by the two large regions of opposite field 
polarity (Fig. 3d), i.e. the field is approximately dipolar. However, 
this structure is noticeably distorted, in particular if one examines 
the distributions of the field modulus and horizontal field (Fig. 3ab). 
Figure 3a suggests that the surface field modulus ranges between 
about 4 and 17 kG. The corresponding disk-integrated mean field 
modulus varies between 6.8 and 7.7 kG depending on the rotational 
phase (somewhat less than was estimated by Glagolevskij & Gerth 
(2008)). 


We also derived a magnetic field map using the inversion 
based on a discrete surface element description of the magnetic 
field. For this inversion to converge to a consistent solution it re¬ 
quired us to initialise the inversions with a basic dipolar field struc¬ 
ture. The resulting fit to the observed Stokes profiles is also shown 
in Figs. 6 and 7. Generally, there is little difference for the Stokes 
IV profiles between these fits and those obtained using the spher¬ 
ical harmonic parameterisation. However, the spherical harmonic 
inversion does a marginally better job at reproducing the Stokes 
QU profiles. A detailed comparison of the radial, meridional and 
azimuthal components of the discrete and spherical harmonic mag¬ 
netic field maps is presented in Fig. 4. It is evident that the inferred 
field topology does not depend significantly on the parameterisa¬ 
tion of the surface field map. 

For an additional test, we derived a magnetic map using only 
Stokes TV profiles. The purpose of this exercise was to assess how 
much of the field complexity is lost in the inversion ignoring lin¬ 
ear polarisation data. The resulting magnetic field map is shown in 
Fig. 4, where it can be directly compared to the results of inversions 
based on all four Stokes parameter spectra. This comparison shows 
that the Stokes IV map does not show as much detail compared to 
the Stokes IQUV maps. This loss of complexity is similar to the 
results of previous MDI studies where the Stokes IQUV and IV 
inversion results were compared (Kochukhov & Wade 2010; Rosen 
et al. 2015). We also find that the meridional field component is 
systematically stronger in the IQUV maps. This is an expected re¬ 
sult as seen from previous numerical tests (Kochukhov & Piskunov 
2002; Rosen & Kochukhov 2012). 

To confirm the reality of smaller scale (£ ^ 3) magnetic fea¬ 
tures in the MDI magnetic map of HD 32633 we investigated if the 
observed Stokes profiles of this star could be reproduced by limit¬ 
ing the code to reconstruct a purely dipole ((Wax = 1 ) or dipole 
+ quadrupole (£ max = 2) field topologies. The best-fitting dipole 
+ quadrupole field distribution is illustrated in Fig. 4. The corre¬ 
sponding fits to observations are shown in Figs. 8 and 9. These fits 
are clearly inferior when compared to those achieved by the inver¬ 
sions which allowed for a more complex field geometry. This result 
demonstrates that the magnetic field of HD 32633 indeed deviates 
significantly from a pure dipole or dipole + quadrupole structure 
and that such low-order multipole models cannot reproduce the ob¬ 
served polarisation spectra. In particular, the Stokes Q and U pro¬ 
files cannot be fit at all. 

The use of the spherical harmonic description of the magnetic 
field enables a detailed characterisation of the field structure at dif¬ 
ferent spatial scales by examining relative magnetic energies of in¬ 
dividual harmonic modes. The magnetic energy distribution for the 
final magnetic map is illustrated in Fig. 5 as a function of the an¬ 
gular degree L The relative harmonic energies also reported in Ta¬ 
ble. 3. We find that the total contribution of the poloidal and toroidal 
field is 84% and 16% respectively. The dipole mode {£ — 1) con¬ 
tains 75% of the magnetic energy; 10% of the energy is in £ = 2 
and the remaining 15% is in £ ^ 3. 


6 CHEMICAL ABUNDANCE DISTRIBUTIONS 

The selection process for determining which chemical elements to 
map was very similar to that used by Silvester et al. (2014b) for 
a 2 CVn. However, unlike in the case of that star, the Stokes I pro¬ 
files of HD 32633 exhibit relatively little variability as a function 
of rotational phase. Therefore lines were selected primarily on the 
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Meridional 


Azimuthal B (kG) 



Longitude (degrees) 


Figure 4. Comparison of the radial, meridional and azimuthal magnetic components for different MDI inversions. From top to bottom: map derived with the 
spherical harmonic inversion and using Stokes IQUV, map derived with the discrete surface element inversion and using Stokes IQUV, map derived with 
the spherical harmonic inversion and using only Stokes IV, map derived assuming a dipole + quadrupole structure and using Stokes IQUV. The dashed line 
indicates the lowest visible latitude for the adopted inclination angle i = 80°. A contour stepping of 2 kG has been used for a range of [—16, +16] kG. 


basis that they did not show significant blending with other atomic 
species. 

The abundance distributions derived here were based on one 
of the following two inversions: abundance maps reconstructed si¬ 
multaneously with the magnetic field geometry (based on spherical 
harmonic parameterisation of the field) and the Stokes IQUV data 
(Si, Cr, Fe and Nd) and abundance maps based only on the Stokes 
TV observations (Mg, Si, Ti and Ni) produced assuming the mag¬ 
netic field structure as derived first from the Stokes IQUV inver¬ 
sion. 

When visualising the final abundance maps in spherical plots, 
the scale was chosen such that extreme outliers were excluded. 
These outliers, typically comprising 2-12% of all surface elements, 
were identified via the use of a histogram of abundance values for 
the entire stellar surface, following the procedure described by Sil¬ 
vester et al. (2014b). 


6.1 Magnesium and Silicon 

The magnesium map was produced using the Stokes IV profiles 
of the Mg II A 4481 feature. The fit between observations and the 
model can be seen in Fig. 11 and the resulting map is presented in 
Fig. 10. The inferred magnesium surface abundance ranges from 
—3.0 dex to —7.0 dex (where the solar value is —4.44) on the 
log(lVMg/iVtot) scale. The high abundance region is located in a 
large belt-like structure above the stellar equator at phases 0.0 to 
0.4 and below the equator at other phases. The high abundance belt 
appears to trace areas where the radial magnetic field is generally 
weak. 

For silicon we derived two maps. One was produced using 
the Stokes IQUV profiles of the Sill A 6347 line and another 
was obtained from the Stokes TV profiles of the Si II lines A 4130 
and 5056. The comparison of observations and the model profiles 
can be seen in Fig. 6, Fig. 7 and Fig. 11. The resulting Si maps 
are shown in Fig. 10. Examining the distribution of silicon in the 
two maps, the abundance values vary from —2.5 dex to —5.5 dex 
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Table 3. Distribution of poloidal (E po \) and toroidal (Et OT ) magnetic field 
energy over different harmonic modes for the best-fitting MDI map of 
HD 32633. 


£ 

Epol 

(%) 

Etor 

(%) 

Etot 

(%) 

1 

69.7 

5.7 

75.4 

2 

9.0 

1.3 

10.3 

3 

1.3 

5.6 

6.9 

4 

1.4 

2.5 

3.9 

5 

0.8 

0.2 

1.0 

6 

0.7 

0.8 

1.5 

7 

0.3 

0.2 

0.4 

8 

0.4 

0.2 

0.6 

9 

0.0 

0.0 

0.1 

10 

0.1 

0.0 

0.1 

Total 

83.5 

16.5 

100 



Figure 5. Distribution of the magnetic energy between different spherical 
harmonic modes in the MDI magnetic held map of HD 32633. Bars of dif¬ 
ferent colour represent the poloidal (blue/dark) and toroidal (red/light) mag¬ 
netic components. 


and from —2.5 to —6.0 dex respectively (where the solar value is 
—4.53). The high abundance areas are located close to the rota¬ 
tional poles, with a large relative underabundance area centred on 
the stellar equator and clearly seen at phase 0.60 in both maps. 

Reconstruction of the silicon abundance distribution from the 
Stokes TV spectra was somewhat complicated by the fact, while 
the A 5056 line could be fit using the YALD oscillator strength, the 
the A 4130 line was under-fit and required the oscillator strength to 
be increased by 0.06 dex. The same problem was encountered by 
Kochukhov et al. (2014) when modelling the same Si II lines in the 
spectrum of Ap star CU Vir. 

Comparing the chemical map with the magnetic field distribu¬ 
tion, in both silicon abundance maps the location of depleted abun¬ 
dance traces areas with weaker magnetic field. On the other hand, 
enhanced abundance areas correspond to regions close to where the 


magnetic field is strongest. At least with regard to the location of ar¬ 
eas of depleted silicon abundance, this behaviour is similar to what 
was found for a 2 CVn (Silvester et al. 2014b). 


6.2 Titanium, Chromium and Iron 

The titanium chemical map was reconstructed using the Stokes TV 
profiles of the Ti II A 4290 and 4563 lines. The fit between obser¬ 
vations and the model can be seen in Fig. 11 and the resulting map 
is presented in Fig. 10. The titanium abundance ranges from —4.0 
dex to —7.0 dex (the solar value is —7.09). The largest areas of 
Ti abundance enhancement are located close to the stellar equator, 
with distinct low abundance regions in the vicinity of the rotational 
poles. The areas of enhanced titanium are located where the field 
strength is somewhat weaker than the average, approximately 4 to 
6 kG. Some of the depleted abundance areas are located close to 
the regions of the strongest magnetic field. 

The chromium map was produced using the Stokes IQUV 
profiles of the Cr II A 4558, 4588 and 4824 lines. The comparison 
between observations and the model profiles is shown in Figs. 6 and 
7. The corresponding map is displayed in Fig. 10. The chromium 
abundance ranges from —3.0 dex to —6.0 dex (solar value is 
—6.40). Similar to what is found for titanium, the high abundance 
regions are located in large structures. In this case, three distinct 
features can be identified, with two located at or close to the stel¬ 
lar equator and one (seen at phase 0.0) extending from the equa¬ 
tor almost to the pole. There are also two small depleted regions 
which are located similarly to the Ti depleted areas. Again, similar 
to Ti, it appears that areas of enhanced chromium are coincident 
with the regions where the field is weaker than the average. The 
small depleted abundance areas are located close to where the field 
is strongest. 

The iron map was obtained from the Stokes IQUV profiles 
of the Fe II A 4583, 4923, 5018 and 5169 lines. The observed and 
synthetic spectra are compared in Figs. 6 and 7, with the resulting 
map shown in Fig. 10. The iron abundance ranges from —2.0 dex 
to —5.0 dex (the solar value is —4.54). Similar to what is found 
for titanium and chromium, the high abundance regions are located 
in large structures close to the stellar equator. Also similar to what 
is seen for titanium, there are smaller depleted regions which are 
located closer to the poles. Comparing the magnetic field structure 
and the iron abundance map, while a clear correlation is absent, 
it does appear that the large enhancement areas correlate with the 
relatively weak field regions while the small depleted abundance 
areas are located close to where the field is strongest. 


6.3 Nickel and Neodymium 

The nickel map was produced using the Stokes TV profiles of the 
Ni II A 4067 line. The observed and best-fitting synthetic profiles 
are presented in Fig. 11 and the resulting map is shown in Fig. 10. 
The nickel abundance ranges from —4.5 dex to —6.5 dex (the solar 
value is —5.82). The map appears to be quite complex. High abun¬ 
dances are localised to small structures, with the biggest structure 
seen close to the stellar equator at phase 0.80. For nickel there is no 
clear correlation between the location of abundance enhancements 
or depletion and the magnetic field structure. 

The neodymium map was derived using the Stokes IQUV 
spectra of Nd III A 5050 and 5102. The fit between observations 
and the model profiles can be seen in Figs. 6 and 7. The resulting 
map is shown in Fig. 10. The neodymium abundance ranges from 
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Figure 7. Same as Fig. 6 for the Stokes QU profiles. 
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Figure 9. Same as Fig. 8 for the Stokes QU profiles. 
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Figure 11. Same as Fig. 6 for the observed (black histogram line) and synthetic (solid blue line) Stokes IV spectra of Mg, Si, Ti and Ni lines used for 
abundance mapping of these elements. 


—5.5 dex to —8.0 dex (solar value is —10.62). The abundance dis¬ 
tribution of this element is fairly homogenous, with enhanced abun¬ 
dance area forming one spot just above the stellar equator at phase 
0.0. This area of Nd enhancement appears to coincide with the re¬ 
gion of large negative radial magnetic field component. 


7 CONCLUSIONS AND DISCUSSION 

In this study we presented the magnetic field maps and chemical 
abundance distributions for the Ap star HD 32633. This is only the 
fourth Ap star for which the surface structure is derived by inter¬ 
preting high-resolution spectropolarimetric data in all four Stokes 
parameters. Our magnetic maps of HD 32633 reveal a largely ax- 
isymmetric, bipolar field topology. The poloidal I — 1 harmonic 
mode is the dominant contributor, with nearly 75% of the energy in 
that mode; 9% of the magnetic energy is contained in the poloidal 
quadrupolar I — 2 mode. We also find 16% of the energy in 
toroidal components. At the same time, we have shown that the 
Stokes parameter profiles of HD 32633 cannot be fit by assuming a 
pure dipole or dipole + quadrupole magnetic field geometry. This 
indicates that the field is definitely more complex than expressed 
by such low-order multipolar parametrisation, despite a small con¬ 
tribution of I ^ 3 modes in the harmonic power spectrum. 

It is useful to compare our MDI results for HD 32633 with the 
field geometries of other A and B stars studied using similar mag¬ 
netic inversion methodologies. In doing so it must be kept in mind 
that MDI limited to Stokes IV data loses some degree of field com¬ 
plexity compared to inversions based on the Stokes IQUV spectra 


(Kochukhov & Wade 2010; Rosen et al. 2015). Therefore, magnetic 
inversions with circular polarisation data provide a reliable infor¬ 
mation about the large scale field components but are less trust¬ 
worthy when it comes to characterising the small-scale (l ^ 3) 
magnetic structures superimposed on the global dipole-like field 
background. 

With these caveats in mind, the late-B star CU Vir 
(HD 124224) mapped using Stokes TV (Kochukhov et al. 2014) 
was found to have 52-60% of the magnetic energy in the I — 1 
harmonic mode, and 23-30% in I — 2 and I — 3. This means that 
the global field topology of CU Vir is somewhat more complex but 
still comparable to that of HD 32633. 

A series of four Stokes parameter MDI studies of a 2 CVn 
(Kochukhov & Wade 2010; Silvester et al. 2014a,b) was performed 
using the version of INVERS10 with the discrete field parameter- 
isation. However, reconstructing the field topology with the input 
data of Silvester et al. (2014b) and the same harmonic code as used 
in the present study, we find 73% of the energy in I — 1 and 10% in 
I — 2, which is almost identical to what was found for HD 32633. 
On the other hand, the total contribution of the toroidal components 
is only ^ 5 % for a 2 CVn compared to 16 % for HD 32633. Thus, 
we find that the degree of complexity of the global poloidal field 
structure of CUVir, a 2 CVn and HD 32633 is similar. All three 
stars have spectral classes B9-A0 and masses in the narrow range 
of 3.0-3.1 M© (Kochukhov & Bagnulo 2006). 

Cooler and less massive Ap stars tend to show simpler mag¬ 
netic field topologies. The field in the extreme Ap star HD 75049 
(Kochukhov et al. 2015) was mapped using Stokes IV and was 
found to be almost entirely poloidal (96% of the magnetic energy 
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in the poloidal harmonic components) and mostly dipolar (90% of 
the energy in £ = 1 mode). The magnetic field of the roAp star 
HD 24712 (DO Eri) was investigated by Rusomarov et al. (2015) 
using Stokes IQUV. They found 96% of the energy in the I — 1 
mode with entirely poloidal structure. Thus, both of these stars have 
a simpler field configuration than found for HD 32633, CU Vir and 
o? CVn. 

Examples of significantly more complex global magnetic 
field configurations are limited to substantially more massive B- 
type stars. For the helium-strong B2 star HD 37776 (V901 Ori) 
mapped using Stokes IV (Kochukhov et al. 2011) it was found 
that the largest contribution to the field topology was in the non- 
axisymmetric £ = 3-4 modes and that the field of that star con¬ 
tained a sizeable toroidal component. The young massive (B0.2V) 
star r Sco (HD 149438), studied with Stokes IV (Donati et al. 
2006), also exhibits a complex field topology, with the £ = 4 
mode being the strongest contributor to the magnetic field structure. 
Another Stokes IV analysis of the B2 He-strong star HD 184927 
(VI671 Cyg) (Yakunin et al. 2015) revealed 38% of the magnetic 
energy in the £ = 1 mode compared to 62% in £ = 2, although 
higher order modes could not be meaningfully characterised for 
this star due to its a low inclination angle. 

A few other A and B-type stars have been mapped with MDI, 
but without the benefit of characterising the field topology with 
spherical harmonics. The magnetic field geometry of 2.1 M© Ap 
star 53 Cam (HD 65339) was studied by Kochukhov et al. (2004) 
using lower quality MuSiCoS four Stokes parameter data. They 
found a dipolar-like field structure distorted by small-scale local 
spots, not unlike the field of a 2 CVn. 

The B2 He-strong star a Ori E (HD 37479) was mapped with 
Stokes IV data by Oksala et al. (2015). This star was found to 
have mostly dipolar magnetic field with a small non-axisymmetric 
quadrupolar contribution. Direct visual comparison of the derived 
magnetic field maps shows that the radial field of a Ori E is less 
complex than that of HD 32633. Overall, the magnetic field geom¬ 
etry of a Ori E is appears to be simpler than the field in HD 37776, 
r Sco and HD 184927. 

The different levels of the global magnetic field complexity of 
A and B stars studied with spherical harmonic MDI are summarised 
in Table 4. The three columns in this table correspond to distinct 
mass ranges. So examining Table 4, it is tempting to conclude, as 
already suggested by Rusomarov et al. (2015), that field complexity 
increased with stellar mass. However, this sample of magnetic stars 
is very small and inhomogeneous in the sense that only three out of 
eight objects were studied in all four Stokes parameters. 

Theoretical modelling by Braithwaite (2009) showed that it is 
possible to have a stable magnetic field configuration for a predom¬ 
inantly poloidal field with some toroidal field components in the 
stellar interior, in essence a bipolar-like magnetic field structures 
as seen in HD 32633 and a 2 CVn. At the same time it has also 
been shown by Braithwaite (2008) that more complex configura¬ 
tions, such as seen in r Sco, could result in a stable field. 

In addition to mapping the magnetic field of HD 32633 we 
also presented chemical abundance maps for 8 elements. Similar 
to what was found in the case of a 2 CVn (Silvester et al. 2014b), 
we find no clear correlation between the location of enhanced and 
depleted abundance areas and the magnetic field. What we do find 
however is that the enhanced areas of iron-peak elements and sili¬ 
con seem to coincide with the areas where the field is weaker, which 
was also found for a 2 CVn. 

Theoretical diffusion modelling of rare earth elements and 
iron-peak elements by Michaud et al. (1981) suggests that rare- 


Table 4. Qualitative comparison of the magnetic field topology complexity 
of HD 32633 with other stars studied with MDI. 


More Complexity 

Similar Complexity 

Less Complexity 

HD 37776 

CUVir 

HD 24712 

T Sco 

a 2 CVn 

HD 75049 

HD 184927 

HD 32633 



earth elements should be concentrated where the magnetic field is 
horizontal and that iron should be enhanced where the field is ver¬ 
tical. Similarly theoretical modelling of silicon by Alecian & Vau- 
clair (1981) expects silicon to be enhanced at the magnetic equator. 
It has been suggested in the case of a magnetic field which was not 
purely dipolar, that instead of belt-like enhancements seen at the 
magnetic equator, only spotlike features should be observed at the 
equator with MDI (Alecian & Stift 2010). None of these theoretical 
predictions are observed in our maps, suggesting further work is re¬ 
quired in the theoretical modelling of diffusion in the atmospheres 
of Ap stars. 
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